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Abstract 


This paper studies the decay of a large, closed domain wall in a closed 
universe. Such walls can form in the presence of a broken, discrete symmetry. 
We introduce a novel process of quantum decay for such a wall, in which the 
vacuum fluctuates from one discrete state to another throughout one half of 
the universe, so that the wall decays into pure field energy. Equivalently, 
the fluctuation can be thought of as the nucleation of a second domain wall 
of zero size, followed by its growth by quantum tunnelling and its collision 
with the first wall, annihilating both. The barrier factor for this quantum 
tunneling is calculated by guessing and verifying a Euclidean instanton for 
the two-wall system. We also discuss the classical origin and evolution of 
closed, topologically spherical domain walls in the early universe, through a 
“budding-off” process involving closed domain walls larger than the Hubble 
radius. This paper is the first of a series on this subject. 
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I. INTRODUCTION 


Following the classic paper of Zel’dovich, Kobsarev & Okun [1]], it lias been generally be¬ 
lieved in early-universe cosmology that domain walls are forbidden, because otherwise they 
would gravitationally dominate the present-day universe, contrary to observation. Conse¬ 
quently, fundamental theories of physics are either forbidden to exhibit spontaneously broken 
discrete symmetries, or are constrained to eliminate the associated domain walls somehow, 
for instance through low energy destabilization of all but one of the discrete vacuua, or 
through decay of the domain walls themselves through nucleation of string loops into walls. 
Only then might an early universe dominated by domain walls evolve into a present-day 
universe consistent with current cosmological observations. 

The purpose of this paper is to discuss a novel process by which domain walls can decay, 
namely, the quantum decay of domain walls by global fluctuation and quantum tunneling. 
We introduce and study this process for closed universes containing a single domain wall, but 
the decay process itself is no doubt more general. Imagine a closed universe separated into 
two different discrete vacuum states by a single closed domain wall. If one of these vacuum 
states spontaneously fluctuates into the other, the domain wall will cease to exist and energy 
will be liberated. Equivalently, the fluctuation can be thought of as the nucleation of a second 
domain wall of zero size, followed by its growth by quantum tunnelling and its collision with 
the first wall, annihilating both. Clearly this decay will be slow since it involves a fluctuation 
over an entire half of the universe, and moreover it involves curved spacetime, and therefore 
requires semiclassical quantum gravity for its description. 

Such a closed universe can be created, for instance, by the the gravitational collapse of 
closed (i.e., bubble-like) domain walls. During a cosmological phase transition that produces 
domain walls by the Kibble mechanism, both open and closed domain walls will be produced. 
Closed domain walls collapse due to their surface tension. A collapsing closed domain 
wall may thermalize completely, leaving no remnant, or it may produce a black hole by 
gravitational collapse. We especially consider in this paper the behavior of closed domain 
walls that happen to be born with size Ro larger than the then-Hubble radius Rh■ These will 
originate naturally with some probability in the Kibble mechanism. As one might expect, 
such a closed domain wall will typically produce a black hole. However, the spacetime 
singularity inside the black hole may not destroy the domain wall; rather, the domain wall 
may expand indefinitely to create a new, inflating universe within the black hole, in a 
“budding-off” process. This process was briefly considered by Blau et al. [§], but has not 
attracted interest since then. Studying a simple but appropriate model, we conclude that 
new-universe creation is obligatory in the model, as long as R 0 > Rh at the time of the 
phase transition, and we argue that this conclusion is apt to be robust. 

The newly created inflating universe is dominated by a single closed domain wall, and 
therefore does not resemble our present universe. Since the domain wall is of finite size, 
however, it is subject, as mentioned above, to decay by quantum tunneling, unlike an infinite 
domain wall, which cannot decay because of an infinite barrier against quantum tunneling. 
After a long time, inflation will therefore end by quantum tunneling of the whole domain 
wall, followed by thermalization of its energy. The resulting universe will therefore exhibit 
a hot big bang, and may conceivably resemble our universe. 

A comprehensive discussion of domain walls in cosmology is given by Vilenkin and Shcl- 
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lard [j3j]. There has been a great deal of work on formation and dynamics of closed domain 
walls between domains of true and false vacuua — the so-called “bubble walls” — includ¬ 
ing the possibility of budding off of new universes within a newly formed black hole. The 
present work differs from that work in that we never appeal to a false vacuum; all domain 
walls in this paper occur between regions of true vacuum, as a result of a spontaneously 
broken, exact, discrete symmetry. There has been some very interesting work on “topo¬ 
logical inflation” within defects including domain walls, cosmic strings or monopoles [f]].^|. 
The present work differs from that in that we consider quantum decay of an entire domain 
wall, while that work depends on quantum fluctuations within a domain wall, monopole, or 
cosmic string. Recently, Caldwell, Chamblin, and Gibbons [6|] studied the pair production 
of black holes in the spacetime of a domain wall, using instanton methods. That process is 
related to, but distinct, from the one we study here. 

In Section II, we discuss some basic properties of domain walls and derive the Euclidean 
action for walls in the thin-wall approximation. Section III reviews the solutions describ¬ 
ing the gravitational fields of open (Vilenkin) and closed (Ipser-Sikivie) domain walls, and 
discusses their behavior. Section IV presents a complete discussion of the gravitational 
fields and motion of domain walls in spherical symmetry, assuming true vacuum everywhere 
except for the domain wall itself. Section V sets up a spherically symmetric cosmological 
model in which a single closed domain wall forms, and demonstrates the conditions under 
which this domain wall can create a new inflating universe. Inflation in this new universe is 
driven by the domain wall that creates it. Section VI contains the heart of our development 
and discusses the quantum decay of the domain wall in a closed, domain wall-dominated 
inflating universe, and estimates the probability per unit time — small but nonvanishing - 
for the decay of the domain wall by tunneling processes in semiclassical quantum gravity. 

We emphasize that all the processes considered in this paper — budding-off, inflation, 
creation of cosmological fluctuations, and quantum wall decay, are driven by domain walls. 
False vacuum energy plays no role. This paper is the first of a series on this subject; the 
second paper will study the same quantum decay process by a different method, namely a 
Hamiltonian approach. 


II. THE ACTION FOR DOMAIN WALLS 


Domain walls appear in matter field theories where discrete symmetries are spontaneously 
broken. The action for matter plus gravity is of the form 


S= d 4 x^/^g 
Jm 


J mat 


+ 


R 


167tG 


+ [ d 3 x ^ Kh . 

JdM 87 tG 


( 1 ) 


where the 4-volume M of the system may have a 3-boundary dM , with 3-metric h and 
3-extrinsic curvature of trace K h ^\ Little of this paper will depend on the details of the 


1 Our conventions are generally those of MTW p|: Greek spacetime indices run over ... = 
0,1,2,3 and the spacetime signature is (—,+,+,+); Roman indices in a hypersurface run over 
a,b,... = 1, 2,3 or a, b ,... = 0,1, 2 for spacelike or timelike hypersurfaces respectively; our extrinsic 
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matter, but for concreteness, the matter may be chosen as a real scalar held (j) with matter 
action 

L mat = ~9 IJV d li (j>d v (l>^U(<j>) (2) 

where U(4>) has a discrete number of degenerate minima at which U = 0, e.g., U((f>) may be 
a “double well” potential. 


A. The Euclidean action 


The Euclidean action / is obtained by analytically continuing Eq. ([[]) to imaginary time 
and then reversing its sign, i.e., 


I 


/ M 


d 4 x^g 


J mat 


+ 


R - 
167rG. 



d 3 x 


VhK h 
8t tG 


( 3 ) 


where the metric is now positive definite. We assume that L mat has a discrete broken 
symmetry and therefore exhibits domain walls as solitons in the low energy limit, without 
any other light degrees of freedom. In this paper, we treat domain walls in the “thin-wall” 
approximation || , where the thickness of the domain wall is assumed much smaller than all 
other length scales in the problem, so that the wall may be taken as a two-dimensional sheet 
of stress-energy in 3-space. The low energy Euclidean action in the thin-wall approximation 
is 


d I'W 


d A x 


v (gR 

167tG 


+ cr 


'Y' [ d 3 xVh — [ d 3 x 

i JDi JdM 


VhK h 

8nG 


( 4 ) 


where the sum on i runs over some finite number of separate domain walls Di. The action of 
each wall is simply proportional to its 3-volume, with the surface energy density a a constant, 
fixed by microphysics of L mat . In the thin-wall approximation, the spacetime geometry is 
generally not smooth across each wall; rather there is a 3-dimensional (5-function in the 
Riemannian curvature at each wall [0. This renders the action (||) awkward to work with. 
Therefore, we will eliminate these (5-functions by breaking up M into a union of 4-dimensional 
voids Mj , each with smooth interior 4-geometry, meeting at domain walls Di. The Euclidean 
action becomes 


t = -£ 


f 

'Mi 167 tG 


y, J d 3 xVh 


a + 


(K-2-K,) 

8nG 


VTiK, 


/ d 3 x 
JdM 8nG 


h 


( 5 ) 


where the prime on the integral f' means that any (5-function in R at a wall does not 
contribute to the 4-integral; and where h a b is now the 3-metric on each domain wall Di (as 
well as the fixed 3-metric on the system boundary dM), and K\ and K 2 are the traces of the 


curvature K a b is the derivative -y£ n h a b of the hypersurface 3-metric h a b with respect to an outward 
pointing normal vector n. Thus, R > 0 and K > 0 for a sphere embedded in flat space. 


4 














extrinsic curvature of each domain wall on its two sides, with respect to the normal vector 
that points from 1 to 2. Due to the unsmoothness of the spacetime geometry, K\ ^ K 2 . 
The instrinsic 3-metric h a b of each domain wall is constrained to agree with that inherited 
from the 4-geometry g /lL/ from the void on each of the two sides. 

The field equations become as follows. Variation of g fiIJ within each M } gives the the 
vacuum Einstein equations 


G^ v = 0 (within Mj) (6) 

within each void. Variation of the domain wall metric h a b on each domain wall Di gives the 
well known Israel jump condition || for the full extrinsic curvature K ab 

Kf - Kf = AnGah ab (7) 

showing that there is a jump 47ra in the trace of the extrinsic curvature at the domain wall, 
while there is no jump in the traceless part of the extrinsic curvature. 


B. The action of a solution 


The solutions of the field equations are the extrema of the action. After the action is 
varied and a solution found, simplifications occur in the formula for value of the action at 
the solution. Substituting the equations of motion (|].[7|) back into the action (|5|) give the 
(weak in the Dirac sense) formula 

hw = [ d^xVh ( 8 ) 

i 2 J Di 


showing that the thin-wall action of a domain wall is always negative. 

This formula was given by Caldwell, Chamblin, and Gibbons [Q for the particular case of 
a real scalar field theory (|]), as follows. The stress-energy tensor for the matter is obtained 
in the canonical way and is found to be 


T = 


2 d(y/^gL mat) 


-<9 M 0<9„0 + g, lb 


-g^d^dpcj) + U (0) 


from which one calculates, taking the trace of the Einstein equations, 

= g a0 d a ^ 30 + AU (0), 

so that the Euclidean action of a solution to general relativity reduces to 


I = -J y/gd*xU(<fi). 


(9) 


( 10 ) 


( 11 ) 


In the thin-wall approximation, U (0) = 0 except very near a wall, where it behaves like a 5- 
function. In a small neighborhood of a point on a wall, the field is approximately planar and 
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is only a function of the coordinate x 3 normal to the wall in a Gaussian normal coordinate 
system. Then the equation of motion for the held has the first integral 


so that Eq. ([|) gives 


1 / d(j)\ 2 

2 


U{<t>) = 0 , 


r; = 2iV)dia g (i, 1,1.0). 


( 12 ) 


(13) 


Furthermore, in this approximation, the stress-energy also has a 5-function singularity at 
the wall [Hip , so that 


T}t = crdiag(l, 1,1, 0)5(x 3 ), (14) 

i.e., 2 U(x 3 ) = a5(x 3 ), where a is the surface energy density of the wall, as well as its tension 
in the tangential directions. Combining this result with Eqs. (10). one again finds Eq.p) 
as the action of a domain wall solution to the held equations. 

For convenience we will adopt the notation 


/I = 4(707 


(15) 


in what follows. 


III. THE VIS SOLUTIONS FOR WALL-DOMINATED UNIVERSES 


A dynamical solution to Einstein’s equations with a single thin domain wall present was 
found by Vilenkin in [IT] and by Ipser and Sikivie in [TT|. We will henceforth refer to this 


solution as the VIS solution. One can construct this solution from hat spacetime 


ds 2 = - dT 2 + dX 2 + dY 2 + dZ 2 (16) 

as follows. The domain wall itself is the hyperboloid 

X 2 + V 2 + Z 2 -T 2 = (J^j , (17) 

and the “gravitational held” of the domain wall is the interior 

X 2 + Y 2 + Z 2 - T 2 < . (18) 

The complete solution is constructed from two such copies of hat spacetime, glued together 
along their hyperboloids. The entire spacetime is therefore spatially closed, with topology S 3 
in the spatial sections. This spacetime can be interpreted as a closed cosmology, containing 
a single spherical domain wall. The universe collapses from infinite size, halts at a minimum 
radius R m i n = due to self-repulsion of the wall, and expands back out to infinite size 
again. 
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The Vilenkin solution is actually a subset of the above S 3 spacetime in a different coor¬ 
dinate system, and moreover it is conventionally given the somewhat different interpretation 
of an infinite planar wall. It is built out of the spacetime: 

ds 2 = dz 2 + z 2 [—dt 2 + exp(fiGt)(dr 2 + r 2 d(j) 2 )\. (19) 

The complete Vilenkin solution is, similarly, two copies of this spacetime for z < glued 
together along the domain wall at z — One can relate the two solutions as follows. 
Rewrite Eq. ( 0 ) slightly as 

ds 2 = -dUdV + dR 2 + R 2 dcj) 2 (20) 


where T = (U + V)/2, Z = (—U + V)/2, X = Rcoscp, and Y = R sin</>. Then the 
transformation relating the solutions is: 


r = 

z = 

t = 



( 21 ) 

( 22 ) 

(23) 


One can verify that the Vilenkin solution given by Eq. (|19|) covers only half of the full 
spacetime, all of which is given by the Ipser and Sikivie result. This is due to the fact 
that the coordinates in (19) cover only half of the (2+l)-dimensional deSitter space in the 
hypersurfaces of constant z. 


IV. MOTION OF SPHERICAL SYMMETRIC DOMAIN WALLS 

In general, domains separated by thin walls may not consist of entirely flat space. In 
this section, we generalize the VIS solution to include nonzero Schwarzschild masses in each 
part of the spacetime, and then analyze the possible resulting dynamics of the wall. 

The junction conditions of the Einstein equations across a thin shell of stress-energy were 
first worked out by Israel in || , who showed that the metric across the surface is continuous, 
whereas the extrinsic curvature has a jump discontinuity there. In the case of spherically 
symmetric domain walls, the trajectory and 4-velocity of the wall are specified by 

x a (r) = (T(t),R(t),6>,0) 

«“(r) = (t(r),R(r), 0,0) 

and the junction condition can be written in the form 

Kjjf 1 - K ( m > = -A GR 2 (26) 

Taking the spacetime on each side of the wall to be Schwarzschild-deSitter, with line element 


(24) 

(25) 
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ds 2 = 
one calculates 


- 1 - 


2 M AR 2 ' 


R 


dT 2 + 1 — 


2 M AR 


2 \ - 1 




df? 2 + R 2 dD 2 , 


(27) 


/ 2M AR 2 ', 
K ee = R 1 - —-—IT 


= ±R 1 


R 3 
2 M AR 2 


R 


+ R 2 


(28) 


where the second equality comes from normalizing the 4-velocity. The junction condition 
can then be rewritten as 


R 2 


R( A+ - A_) (Af+ - M_)l 2 (/iGR) 2 

6/1 + ^R 2 J + 4 

G(M + + M_) GR 2 (A + + A_) 

R + 6 ’ 


(29) 


which is a first order ODE for R(t), and can be regarded as an energy-type equation with 
p, A and M all conserved. The second order equation of motion for R(r) can be found by 
differentiation. Restricting to the true vacuum, one finds 


R 2 = 


(M+ - M_) 

/IR2 


(/iGR) 2 G(M + + M_) 
+ 4 + R 


1, 


(30) 


and further restricting to the case of flat spacetime, one finds 


R 2 


(nGR ) 2 

4 


(31) 


which is the VIS solution. In what follows we will take A = 0, and study the behavior of 
domain walls in the absence of vacuum energy, i.e., with dynamics given by Eq. (|30|). 

To examine Eq. (j3(]) qualitatively, change to the dimensionless variables 


z = R 


H 2 G 


1 1/3 


r = 


2(M+ - M_) 

fiGr 


upon which the equation of motion of the wall becomes 

z' 2 + V(z) = E, 


(32) 

(33) 


(34) 


where 


V{z) 

E 


z 6 + 2 z 3 (M + + M_)/(M+ - M_) + 1 


/liG 2 (M+ - M_) 


(35) 

(36) 
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Although we will assume M + > M _ so that z > 0 in what follows, this is a matter only of 
convention since Eq. (30) is symmetric with respect to M + and M_. 

Figure 1 shows a typical plot of V(z), and one sees that there exist several qualitatively 
distinct classes of solutions, depending upon the values of the parameters (M_, M + , ft). A 
wall may be born with zero size, expand to a finite maximum radius and recollapse; it may be 
born with zero size and expand indefinitely; it may collapse from infinite size to a minimum 
radius and then reexpand (this is the behavior of the VIS solution); or it may collapse from 
infinite size to zero size. In each case, the complete spacetime consists of two pieces of the 
extended Schwarzschild solution, glued together at the wall. 


9 



FIGURES 



FIG. 1. The “potential energy” function V{z) for domain walls moving in Schwarzschild space- 
time. When a trajectory crosses the vertical dashed line at z = 1, then the polar angle (angle 
between the origin and the trajectory) in one half of the complete spacetime changes sign. The 
horizontal lines represent the possible qualitatively distinct classes of dynamics. 


One can construct these spacetimes as follows. Change to coordinates which are well- 
behaved at the Schwarzschild horizon, and in which light rays travel at 45 degrees. Then 
there is a correspondence between the sign of the extrinsic curvature in Eq. (|28|) and the sign 
of the wall trajectory’s “angular velocity”, ^ arctan(a: 0 /a: 1 ) in those coordinates. For exam¬ 
ple, in the Kruskal-Szekeres B coordinates ( v,u,9,(j )), which are related to the standard 
Schwarzschild coordinates by 


one finds 


u 


T 



R 

g 2 M 


4Marctanh 
4Marctanh (j±\ 


R > 2M, 
R < 2M, 


K+ e = ±R 




Kee = 


+8 M\e 

±R (l - 


R 

2M + (uv — vii ); 


2 M_ 
R 



1 

2 


(37) 

(38) 


(39) 
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= —8M^e 2M - (mu — ra). (40) 

In these coordinates, the polar angle in the spacetime diagram is given by tan Q — v/u, and 
its rate of change is given by 


c 


uv — vu 
u 2 + v 2 


(41) 


Hence we see that when K^ g is positive in a given region of the spacetime, the wall trajectory 
follows a path such that Q increases along the trajectory, and vice versa. Similarly one sees 
that when K go is positive, ( decreases along the wall trajectory, and vice versa. The sign 
change is due to the fact that the vector normal to the 4-velocity changes its orientation 
from one side of the wall to the other. 

The last step is to determine the signs of the quantities K^ e and K gg as a function of z 
in Fig. 1. One finds that 


whereas 


Kte 


> 0 , 

< 0 , 


z > 1 
z < 1 


Kqq <1 0, z 0. 


(42) 


(43) 


Finally one knows everything necessary to construct the complete spacetimes from Fig. 1. 
In each part of the spacetime, one only has to draw the wall trajectory such that it has the 
correct starting and ending points, and also so that the angle ( increases or decreases in 
each part according to the signs of the extrinsic curvature as stated above. 

Figure 2 shows the Penrose (conformal) spacetime diagrams for each of the possible 
cases. In each case, the complete spacetime consists of the wall trajectory, those points in 
the right-hand copy of Schwarzschild to its right, and those points in the left-hand copy 
of Schwarzschild to its left. Of particular interest is case (iv), where a wall is born with 
zero size and grows to infinite size. An observer in the asymptotically flat region IV of 
the spacetime would see the wall enter the Schwarzschild horizon; however, the wall subse¬ 
quently avoids the singularity within the horizon and expands indefinitely, creating a new, 
inflating universe inside the black hole. Clearly, this model suffers from the requirement 
that an initial singularity exist from which the domain wall emerges. We therefore proceed 
to explicitly construct a model which circumvents this necessity, while still producing a new 
wall-dominated universe. 
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V. 




FIG. 2. Spacetime diagrams for each of the possible qualitatively distinct classes of dynamics, 
for domain walls moving in Schwarzschild spacetime. In each case, the complete spacetime consists 
of the wall, those points in the right-hand copy of Schwarzschild to its right, and those points in 
the left-hand copy of Schwarzschild to its left. 


V. A MODEL OF NEWLY FORMED BUBBLES OF DOMAIN WALL 

We now set up and study a toy model of the behavior of topologically closed and spher¬ 
ical domain walls in cosmology. For convenience we work in units where G — 1, in this 
Section only. Our cosmological model is spherically symmetrical, closed, and bounded. The 
matter content of the model is pressureless dust, arranged in concentric spherical shells. An 
important feature of our model is that would be fated to recollapse in a finite time, if it 
were not for the newly born closed domain wall. This feature is important because we want 
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to avoid giving the newborn domain wall a spurious kick, i.e., a positive total energy. In 
applications we may think of the model as being only moderately larger than the Planck 
radius, and with a total lifetime only moderately larger than the Planck time — if it were 
not for domain wall effects. 

(As a limiting case we can study a universe in which each shell starts with marginal 
(“/c = 0”) binding energy, so that the whole universe just barely expands to infinity as 
t —> oo.) 

In other words, we want our model to be a purely curvature-free pertubation of a com¬ 
pletely homogeneous pressureless k = +1 model, because phase transitions cannot produce 
correlated pertubations outside the light cone, and therefore, as is well known, cannot pro¬ 
duce curvature perturbations. 

We concentrate on a single dust shell of finite total proper mass AM, idealized as being 
infinitely thin. J] Spacetime is vacuum for a small region (see below) around our shell, to 
avoid complications due to collisions of shells; see below. By Birkhoff’s theorem, the space- 
time just inside the shell is Schwarzschild of some mass M l5 and that outside is Schwarzschild 
of some mass M 2 . Our shell moves according to the first-order equation of motion 


AM 2 M 2 + Mi 
4R 2 + R + AM 2 


(44) 


where i?(r) is the curvature radius of our shell as a function of its proper time r, and =d/dr. 
In order that the shell be gravitationally bound, with some relative binding energy e, we 
take 


M 2 - Mi = eAM (0 < e < 1). 


Then the equation of motion becomes 


R 2 


AM 2 M 2 + Mi 
4 R 2 + R 


(45) 


(46) 


If 0 < e < 1 our shell would recollapse; in the limit that e = 1 our shell would just make it 
to infinity — if left undisturbed. However, it is instead going to turn into a domain wall. 

To model the birth of a domain wall, our shell is assumed to suddenly become a domain 
wall as it reaches some radius Rq: a tangental stress equal in magnitude to the surface 
mass density suddenly appears. However, the mass density (and total mass AM) must 
remain meanwhile remain constant, according to conservation laws. Thereafter, our shcll- 
cum-domain-wall moves according to a different first-order equation of motion, 


r.2 1.2 02 , , (M2 + M,) , (M 2 - 

R =-^R -1 +- 5 


(47) 


2 An infinitely thin “dust” shell of finite mass is an acceptable idealization within general relativity; 
however, subtleties do lurk: A close analysis shows that the shell must possess radial stresses in 
order to stay infinitely thin — this is for instance clear from from the fact that a dust particle just 
inside the shell must experience a different acceleration than one just outside. 
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( 48 ) 


To match R and R at Rq (these must be continuous) we have 

AM 


fi = 


«0 


Thus the first-order equation of motion of R for the spherical domain wall becomes 


R 2 


H 2 R 2 M 2 + M 1 e 2 Rt 

4 + R + i? 4 


(49) 


The main question is simple: How does the solution of this equation R(r) behave? Does 
it go to ex) as r —> oo, or in contrast does it fall back to R = 0 at some finite time r? 
The former behavior means that the shell has formed a new, inflating, semi-closed universe 
inside the black hole M 2 , while the latter behavior means that the shell has crunched into 
the singularity inside the black hole M 2 . Note that in both cases, the domain wall creates 
a new black hole of mass M 2 by gravitational collapse, as seen by external observers in the 
original universe. 

Our model therefore has three dimensionful paramaters, /i, M \, and Rq, and one further 
dimensionless parameter, e. (By the above, M 2 = M\ + neR 2 .) Physically, // is fixed by 
microphysics, Ro is set through the Kibble mechanism by random variations during the 
phase transition and is likely to be roughly the Hubble radius then, and M\ is the mass of 
matter that happens to lie inside the bubble of domain wall when it forms. The relative 
binding energy e is set by the overall dynamics of the initial universe, with some further 
adjustment by the Kibble mechanism. 

Since three of the parameters are dimensionful, it is convenient to render them dimen¬ 
sionless by forming ratios with the quantity M 2 + Mi, reducing them to two dimensionless 
parameters, m and r 0 , with 0 < m < 1 and 0 < r 0 < 00 : 


M 2 — M 1 fit /? ( j 

772 — -—- = - 

M 2 + Mi jieRl + 2 Mi ’ 

_ i?0 _ -Ro 

r ° ““ M 2 + Mi ~ /ieRq + 2 Mi ’ 

and a dimensionless radius r(t) as a function of a dimensionless proper time t, 

R _ R 

r ~ M 2 + Mi ~ neRl + 2Mi ’ 
r r 

~ M 2 + Mi ~ iieKl + 2Mi ‘ 


(50) 

(51) 

(52) 

(53) 


In terms of the dimensionless quantities, the equation of motion is 

r' 2 + V (r) = 0 


where the effective potential is 


V (r) 


m 2 o 1 

- r 2 + 1 — — 

4 e 2, r 4 r 


2,„4 


6T, 


(54) 


(55) 


The effective potential V(r) diverges to —00 as either r—> 0 or r —> 00 . Therefore V(r) has 
a maximum at some point r x which will be studied below; it will turn out that V(r) has 
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only one local maximum for 0 < r < oo. If the system point gets to sufficiently large r, 
a new, semi-closed, inflating universe will be born. The system point starts at r = r 0 , at 
the moment when the domain wall is born. The key question is whether it can get past the 
maximum at r x . It will get to large r if either of two conditions is satisfied: 1) r 0 > r x ; or 
2) V(r x ) < 0. 

Let us find the maximum of V(r). The gradient of V is 

dv/dr = -2&? + 1 /r2 + ^r (56) 

and the only positive root of the equation dV/dr = 0 is found to be 


4/3 2/3 _ wo 

r - = -^7 w ( x + Vl + 8m2 ) • 

The condition 1) above that r 0 > r x is then 

9 

m 

° — e 2 (l + \/l + 8m 2 ) 

The maximum value of the potential is 

3 m 2 / 3 1 + 2 m 2 + \/l + 8 m 2 
2e 2 / 3 rt /3 (l + Vl + 8m 2 ) 4/3 

and condition 2) that this maximum be negative is 


r o < 


3 \ 3/4 m 1/2 (l + 2m 2 + y/l + 8m 2 ) / 
2/ e 1 / 2 1 + a/1 + 8-m 2 


(57) 


(58) 


(59) 


(60) 


Only one of Eqs. (|58| , [6^ ) need be satished for the system to get to large r. Comparing these 
equations one sees that Eq. (^j) dominates (he., is less restrictive than) Eq. (|58p as long as 


e > 


2^ 1//2 m 

3/ (1 + 2 m 2 + V1 + 8n i 2 ) 1 / 2 


(61) 


while on the other hand Eq. (|58|) dominates if this inequality is reversed. In particular, Eq. 
(1611) holds for all permissible m if e > 1/3. So let us concentrate on the case e > 1/3 (shell 
not greatly bound) and study Eq. (6fj). 

Returning to dimensional variables, the condition Eq. (|60|) becomes 


e/i 2 i?o + 2/iMi > 


2\ 3 / 2 (\/l + 8m 2 ) 

^ ^1 + 2m 2 + \/l + 8m 2 ) ^ 


(62) 


This equation is the condition that the newly formed bubble will form an inflating, semi- 
closed universe. The RHS of Eq. ( |62|) is of order unity and nearly constant for 0 < m < 1, 
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decreasing from 0.77 at m = 0 to 0.59 at m = 1. In particular, this equation shows that 
such a universe will always form if the initial bubble radius Rq is large enough, 

fiR 0 > 0.8/e 1 / 2 . (63) 

Furthermore, Eq. (|62| ) can be shown to imply that the bubble is outside the Hubble radius 
when it forms, in the sense that ‘IM^/R < 1. 

Equation (|62| ) is also sufficient to form a new inflating universe even if e is small and Eq. 
(|6l|) is not satisfied, but it is not necessary. In that case it is easier to form a new inflating 
universe than implied by Eq. (|6^). 


VI. INSTANTONS AND QUANTUM DECAY 


We now wish to estimate the probability per unit time that the wall-dominated universe 
we have described in the previous section will decay due to quantum effects. For simplicity, 
we will consider the case of the VIS solution; however, in the absence of false vacuum energy, 
one would expect the effect of a nonvanishing Schwarzschild mass to be small. We proceed 
by constructing a Euclidean instanton which can be interpreted as interpolating between a 
time slice of the VIS solution and a time slice of a two-wall spacetime in which a second 
domain wall has nucleated at zero size, with the two walls tunnelling toward each other until 
they collide and annihilate. 

We begin by calculating the action of the Euclidean VIS solution. Recall that the radius 
of the wall as a function of the Minkowski time in flat coordinates is given by Eq. m 

X 2 + Y 2 + Z 2 -T 2 = , (64) 


and that the complete spacetime consists of the interiors of two such hyperboloids, identified 
at their boundaries (the wall). The Euclidean VIS solution therefore consists of 2 balls of 
flat 4-space, with the wall as their common (identified) boundary. Then by Eq. (|8|), the 
action of this solution is 


Ivis = J Vhd 3 x 

= -|( 2 ") 

2t r 


(65) 


Having constructed the Euclidean VIS thin-wall solution, we can now construct a more 
general class of thin-wall Euclidean spacetimes, some members of which are solutions to the 
held equations, and therefore candidates for instantons that mediate interesting processes 
such as quantum tunneling. We construct these spacetimes out of some even number 2 n 
(n = 1, 2, 3...) of pieces; the Euclidean VIS solution, composed of two pieces, will be the 
case n — 1. Each piece is a lens-shaped region of hat 4-space bounded by two 3-spherical 
segments. The two 3-spherical segments join on a complete 2-sphere (“the edge of the lens”). 
We form the complete space by identifying pairwise the 3-spherical segments in round-robin 
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fashion; all the 2-spherical edges are identified to a single 2-sphere. We demand that space 
be locally flat at this 2-sphere (no conical singularity — angle of 2n around any circle). In 
the whole space, the walls are precisely the 2 n identified 3-spherical segments. It will now 
be shown that for a given value of n, a configuration which extremizes the Euclidean action 
consists of 2 n identical lenses, with the two 3-spherical segments bounding each lens meeting 
at an angle 27i/2n. Furthermore, the radius of curvature of each segment will be equal to 
2/fiG, the radius of the VIS solution. 

For a flat, compact spacetime, the action in the form of Eq. (||) becomes 

hw = f d 3 xVh 

i J Di 


~{K i + A 2 ) + 2/i 

87 tG 


( 66 ) 


where the sum is over the 2 n 3-spherical lens boundaries as described above. Each 3-spherical 
segment can be described as a segment, or cap, of S 3 of some fixed radius a* and maximum 
polar angle Of 

ds 2 = a 2 [d9 2 + sin 2 OdXl 2 ^ , (67) 

where 0 < 0 < 0 { . 

The extrinsic curvature and its trace at the surface of a 3-sphere of radius a are given by 
K 1 1 = K 2 2 = K 33 = ±l/a and K = ±3 /a, where the sign depends on whether the sphere 
is chosen to have positive or negative curvature on each side. The case of a positive-energy 
domain wall corresponds to choosing negative curvature on each side of the wall, i.e., an 
observer on either side of the spherical wall is enclosed by it. Using K\ + K 2 = — 6 /a in the 
above action, we arrive at 


hw 4tt?. 


d 3 xVh 


Di 


-3 


_G&i 


+ V* 


( 68 ) 


In order that the geometry be smooth, we demand that the area of the 2-sphere which joins 
any pair of 3-sperical segments be the same for all pairs of segments, hence 

A 2 = J V9xx9H d xd^ = 47 ra 2 sin 2 6 *j = const = 47 ra 2 , (69) 

where a is a constant for all lenses in a configuration. Likewise we demand smoothness (no 
conical singularity) at the pole of our configuration; hence 


2 n 

= (TO) 

i=l 

These constraints can be used to eliminate the 2 n radii a* in favor of the single parameter 
a, and to eliminate one of the polar angles in favor of the remaining 2n — 1. The physical 
parameter space is now (0 < < 7r, a > 0), where i — 1, ..., 2n — 1. 

The action is then 


7n(di, • • •, 0 2n — 1 , a) 

1 2 ^ 7 x / sin 26i \ / /ia 3 3a 2 \ 

2 V ' 2 ) y sin 3 di ~ G sin 2 oj 

l( _ , sin2E^A / Ata 3 3a 2 \ 

+ 2l v 7r ' 2 J Un 3 £^ Gsin 2 E^J‘ 


( 71 ) 
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To extremize the action over our configurations, we set 


0 = 


9I n -2 a 

Sft = sm 


n 


f a \ 3 3 fa 

Vsin Oj 


3 cos 9i ( sin 26 i ' 

6i --— 


G V sin 9 t 
2 / a N 2 


2 sin 9i 
— sin 2 Ed 




a 


G Vsin 9i 
3 / a 


-li 


sin 9i 


. sin E 9i > 

2n—l 


3 cos E9i f ^ 
2 sm E Oi V 


G V sin E 9% 
sin 2E0<\ 


(72) 


X 


G Vsin E 9i 


-/i 


E sin 9i 


and 


0 = 


d I n 

da 


2n—l 




'Oi 


| sin 29i ' 


2 ^ \ sin 9, 
1=1 \ 1 


3 / 7r-E^i + ^sin2E^i 


sin E Oi 



(73) 


which has the obvious, symmetrical solution 


9 1 = .. 

a 


sin 9j 


— #2n—1 
2 

JiG' 


7T 

2n’ 


(74) 

(75) 


We believe that there are no other allowed solutions, but have not proved this. Although 
we have extremized the action only within our set of configurations, in fact the extremizing 
spacetimes are solutions of the full Euclidean field equations. Next we determine whether 
these extrema are minima, maxima or saddlepoints of the action. 

To determine the character of an extremum of a function of several variables, one con¬ 
structs the Hessian matrix, or matrix of second partial derivatives. Then the extremum is a 
local minimum if and only if every eigenvalue of that matrix is nonnegative, whereas each 
negative eigenvalue represents a direction in which the function decreases, so that if such 
eigenvalues exist then the extremum is a saddlepoint [jnj. In the present case one finds 


Hd^ —Ax 

Haa = B, 

He ia = 0, 


f~ 2 
-1 


-1 

-2 


V-l -1 


... -IX 
... -1 

-1 - 2 / 


(76) 

(77) 

(78) 


where A and B are positive constants for a given value of n, defined by 
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A = 


2 

p 2 G 3 


„ . 7T „ 9 7T 

( 7T 

. 7T \ 

2 sm-3 cot 2 — 


— Sill — 

n 2 n 

Vrt 

nj 



7 r — n sin - 

n 


sin 2 


7T 

2n 


(79) 

(80) 


The angular part of the Hessian has eigenvalues (— 2n, —1,..., —1), so that this extremum 
is clearly a saddlepoint of the Euclidean action. However, since no other extrema appear 
to exist in the physically allowed parameter space, this symmetric solution can still be 
interpreted as an instanton which mediates quantum tunneling.]] 

The total Euclidean action of this solution is 


In 


2i r 

p 2 G 3 



(81) 


Since this action is an increasing function of n, the VIS solution (n = 1) is the thin-wall 
solution of least action, and the next lowest lying solution will be n = 2, which we therefore 
choose as our candidate instanton. It consists of 4 lenses of flat 4-space, bounded by segments 
of 3-spheres which we identify pairwise and each pair of which meets at the angle 7 t/ 2. One 
can think of this instanton as mediating the creation from nothing of two domain walls; 
from Eq. (]69|), each wall has radius of curvature a 2 = cq sin(7r/4), where aq = 2/(/iG) is the 
VIS radius. Hence the two walls are created by the instanton with radius R = \/2 /(pG). 

One can choose a “final” slice of zero extrinsic curvature through this 4-geometry such 
that at an instant of Euclidean time, a single wall separates two domains which contain the 
same phase of the scalar field; the slice passes through the centers of two of the lenses as 
shown in Figure 3. Let us refer to these two lenses as the primary lenses. Then the slices 
evolve backwards in Euclidean time towards the VIS solution as follows: successive slices 
intersect the primary lenses in a sequence of 2-spheres of decreasing radii, with each slice 
passing through one of the intermediate lenses in such a way that the 2-spheres shrink faster 
in one of the primary lenses than in the other. This process continues until the smaller 
of the 2-spherical slices has reached zero radius, as shown in Figure 4; this is the point at 
which the second wall appears. We require that the initial slice be isomorphic to one of zero 
extrinsic curvature through the VIS space; however, our instanton contains no such slice. 


3 In general, it is not possible to determine absolutely that stable solutions to the nonlinear 
Eqs. ( [VID do not exist in the physical parameter space (0 < di < it, a > 0), i.e., that a local 
minimum of the action does not exist and make the dominant contribution to the Euclidean path 
integral for the tunneling amplitude. However, we have searched for such solutions numerically, 
and we do not believe that they exist. 
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FIG. 3. A hypersurface of zero extrinsic curvature passing through the n = 2 symmetric in- 
stanton. The surface contains two regions of space where the phase of the scalar field is the same. 


8 8 





FIG. 4. The final stages of the evolution of the hypersurface which interpolates between the final 
and initial configurations. Successive instants of Euclidean time are numbered; slice 5 represents 
the instant at which the second wall is created with zero radius. 

Although the meaning of this is unclear, similar pathologies have been encountered in 
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Euclidean quantum gravity before. One way of getting around the difficulty is simply to 
adopt the no-boundary proposal, i.e., to calculate the barrier factors for tunneling from 
nothing of the initial and final held configurations, and then take their difference as the 
action of the interpolation. This was the approach taken, for example, in Hl4|,[T5|,|6| . Another 
approach has been proposed by Farlii, Guth and Guven [JT6[] , who encountered this difficulty 
when calculating the amplitude for a false-vacuum bubble to tunnel through a classically 
forbidden region, so that it could expand perpetually. They suggested the following rule: 
that the Euclidean manifold exhibits a 2-sheeted structure, with the interpolating slices 
moving for part of their evolution on the second sheet, which would then contain a suitable 
initial slice. (One must then fold a single slice which can be matched between the two sheets.) 
In their work, they called such a multi-sheeted Euclidean manifold a “pseudomanifold”, Q 
and they defined the covering number of any point in the manifold as the number of times 
the point is crossed by the evolving hypersurface in the future direction, minus the number 
of times it is crossed in the past direction. They then found that the action weighted by the 
covering number yields the correct Euclidean equations of motion. 

In the present case the second sheet, on which the slices evolve before the new domain 
wall has appeared, is simply the Euclidean VIS manifold. The evolution then proceeds as 
shown in Figure 5. The action of the interpolation is calculated by following the hypersurface 
through its complete evolution, and summing the wall area (and thus the action) “swept 
out” by it, weighted by the covering number as described above. For the first part of the 
evolution (slices 1-4 in Fig. 5), the action is 


-Ivis - A/, 


(82) 


where Iyis is the Euclidean action of the VIS solution, and A/ is the action of the area of 
the 3-hemisphere not swept out by this part of the evolution. Now the slice begins evolving 
on that part of the manifold containing the instanton, and the action of this part of the 
interpolation is 



(83) 


where Jo is the Euclidean action of the instanton, and the minus sign is present because 
the slice evolves through these points with the opposite orientation. Summing the two 
contributions, we see that the total Euclidean action of the interpolation is 

I = 2^vis ~ h)- (84) 


4 Not to be confused with a pseudomanifold in topology. 
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FIG. 5. The initial stages of the evolution of the interpolating hypersurface, which take place 
on the VIS sheet of the manifold. The second sheet, consisting of the n = 2 instanton, is indicated 
by the curved, dashed lines. 


A physical interpretation of this process is that in one half of the VIS space, a second 
wall is nucleated quantum mechanically and grows in imaginary time while the first wall 
shrinks; when the walls reach the same place and same size, they annihilate in a burst of 
energy. In other words, one imagines the scalar held in one half of the space fluctuating, in 
its entirety, into the opposite vacuum state. At the same time, the geometry fluctuates so 
that the domain wall radius is reduced by a factor of y/2 at the instant of its decay. 

The Euclidean action of the instanton is 



4 - 2vr 


/i 2 G 3 ' 

The action of the interpolation is therefore 



I 


1 

2 


{Wis ~ h) 


2 

/i 2 G r 


(85) 


( 86 ) 


The amplitude for the decay process is then given in the semiclassical limit by |L7 


T 

V 


Aexp[—/], 


(87) 
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where the prefactor A can by calculated by considering perturbations about the instanton 
0. Ignoring the prefactor, we thus find that the probability for decay will contain the 
barrier factor 


P ~ exp 


2 

/i 2 G 3 


( 88 ) 


VII. CONCLUSIONS 


The study of domain walls in the early universe is motivated both by their inflationary 
nature, and by the naturalness with which they arise from fundamental held theories. The 
spacetimes associated with domain walls are in many ways analagous to deSitter spacetime, 
with the vacuum energy confined to the two-dimensional sheets which are the walls. How¬ 
ever, the presence of domain walls does not require a global source of vacuum energy; only 
a broken discrete symmetry need exist to lead to the formation of the walls. Domain walls 
may arise more naturally in microphysics than the false vacuum which commonly drives 
models of inflation. 

As in the case of false vacuum energy, domain walls tend to dominate the classical dynam¬ 
ics of the universe due to their gravitational properties. However, the resulting cosmological 
models exhibit large-scale anisotropies in the CMBR, inconsistent with observations, and 
although the instability of the false vacuum is well-known [|17[] , there was until now no corre¬ 
sponding result in the case of pure domain walls. For this reason, it has long been believed 
that domain walls should be forbidden in the early universe. The primary result of this 
paper has been to demonstrate that closed domain walls are in fact quantum mechanically 
unstable, and will therefore decay with finite probability. 

Here we have focused on the evolution of closed, spherically symmetric domain walls, 
which in general have a Schwarzschild mass and an associated singularity. In our analysis 
of the classical dynamics of these spacetimes, we have shown that four general classes of 
behavior are possible: a wall may be born with zero size, expand to a finite maximum 
radius and recollapse; it may be born with zero size and expand indefinitely; it may collapse 
from infinite size to a minimum radius and then reexpand; or it may collapse from infinite 
size to zero size. In the case where the domain wall expands indefinitely, it passes inside 
the black hole horizon, avoids the singularity and creates a new, inflating universe which 
is causally disconnected from the original spacetime, as first pointed out by Blau et al. [Q. 
We have shown that this process will occur naturally if the spherical domain wall is formed 
with radius larger than the Hubble radius at formation. 

Our toy model to study the decay of such a domain wall-dominated universe was chosen 
to be the limiting case where the Schwarzschild mass of the wall vanishes. In this case, we 
have seen that there is an instanton which mediates the decay process, and that the decay 
probability per unit time is 


P 


exp 


/i 2 G 3 


(89) 


Although the process is heavily suppressed, we emphasize that in our scenario there are no 
competing processes; the domain wall simply expands until it decays, whenever that may 
be. The energy of the wall will subsequently become thermalized and lead to a hot big bang. 
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Our instanton calculation also predicts the classical state immediately after the decay, 
namely a closed universe of smaller volume, consisting of two regions of flat space meeting 
at a wall of energy. The wall of energy is the annihilation product of the two domain walls; 
its geometry is that of a 2-sphere with a radius y/2//iG, or a factor l/\/2 smaller than the 
radius of the domain wall before the decay. 

Discussion of the subsequent fate of this new universe is beyond the scope of this paper. 
But this closed universe now contains a form of energy which obeys all the usual energy 
conditions, and therefore recollapse may be expected. If the wall of energy is in a metastable 
state (“slow rolldown”) then recollapse may be delayed greatly, however, and this universe 
might expand to a much greater radius before recollapse. 

Our instanton calculation has also raised some interesting questions about the overall 
nature and validity of the instanton approach to quantum gravity. In general, one would 
like to find a Euclidean instanton which, in addition to being the least-action solution of the 
Euclidean held equations, satisfies two conditions: (i) it contains slices which are isomorphic 
to static slices in both the final and the initial state of the tunneling process, and (ii) the 
hypersurface which interpolates between these two slices has a unique trajectory, ffowever, 
as in our calculation, one Ends that this is not always possible; in such cases, it is not clear 
how to proceed. We have followed the rule of Far hi, Guth and Guven [fid) 
interpolating hypersurface on a 2-sheeted 


to evolve the 


“pseudomanifold” made of two instantons glued 
together over a common region, wherein the Erst sheet contains the final slice, and the 
second sheet which contains the initial slice. The two parts of the evolution are joined 
on a surface in the common region. An alternate approach, which has become standard 
practice (see, for example, HI), is simply to subtract the action of the instanton which 
mediates “tunneling from nothing” of the initial state, from that of the instanton which 
describes “tunneling from nothing” of the final state, and take this result as the action of 
the interpolation. This approach does not even require that the two instantons have a single 
hypersurface which can be identified. The two approaches give identical answers in the 
present case, although the reasons for this are unclear. 

Furthermore, our instanton is a saddlepoint, rather than a minimum, of the Euclidean 
action for the class of field configurations we have considered, as invariably happens in 
quantum gravity. Since there do not appear to be any extrema of lesser action, we have 
taken our solution to be the dominant contribution to the decay process, ffowever, we have 
not proven that this is the case. 

In order to avoid the ambiguities in the instanton calculation, we turn directly to the 
canonical quantization of the the minisuperspace model corresponding to the domain wall 
spacetime under consideration. Since the physical interpretation of the decay process is that 
of a second domain wall being nucleated in one half of the existing spacetime, and the two 
walls tunneling towards each other, there will be two degrees of freedom in the model. We 


will carry out the canonical quantization of such a two-wall spacetime in the next paper [JT8 
and we will find a significantly different result for the quantum decay. 
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